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We consider an extended model of DBI massive gravity by generalizing the fiducial metric to 
be an induced metric on the brane corresponding to a domain wall moving in five-dimensional 
Schwarzschild-Anti-de Sitter spacetime. The model admits all solutions of FLRW metric including 
flat, closed and open geometries while the original one does not. The background solutions can 
be divided into two branches namely self-accelerating branch and normal branch. For the self- 
accelerating branch, the graviton mass plays the role of cosmological constant to drive the late-time 
acceleration of the universe. It is found that the number degrees of freedom of gravitational sector 
is not correct similar to the original DBI massive gravity. There are only two propagating degrees of 
freedom from tensor modes. For normal branch, we restrict our attention to a particular class of the 
solutions which provides an accelerated expansion of the universe. It is found that the number of 
degrees of freedom in the model is correct. However, at least one of them is ghost degree of freedom 
which always present at small scale implying that the theory is not stable. 


I. INTRODUCTION 

Massive gravity is a theory of a massive spin-2 graviton, a generalization from the Einstein gravity which corresponds 
to the massless graviton. The linearized massive gravity was constructed by Fierz and Pauli [1] by adding an interaction 
term into the linearlized Einstein-Hilbert action. However, the Fierz-Pauli theory suffers from the disagreement 
between the predictions made from the massless limit of the theory and those made from the general relativity, 
which was proposed by van Dam, Veltmam and Zakharov (known as vDVZ discontinuity) Si- The solution to the 
discontinuity was clarified by Vainshtein [^, that the linear approximation breaks down for the massless limit of the 
theory and then such discontinuity can be lifted in non-linear theory of massive gravity. The cost of the non-linear 
generalization is the existence of the ghostly sixth degree of freedom, the Boulware-Deser (BD) ghost [1|. The ghost 
free generalization of massive gravity, the so-called dRGT massive gravity, was successfully constructed by de Rham, 
Gabadadze, and Tolley [1,0 which the suitable interaction terms are chosen such that there exists no BD ghost in 
the theory. 

In order to investigate the cosmological implication of the dRGT massive gravity, it is convenient to use the 
Friedman-Lemaitre-Robertson-Walker (FLRW) metric as the physical metric. It was found that the dRGT massive 
gravity does not admit flat and closed FLRW solution iS while the self-accelerating open FLRW solutions were 
found in 0. A more general issue on the model was investigated by replacing the Minkowski fiducial metric with a 
de Sitter or FLRW fiducial metric [Tol - [H |. It was found that the model admits not only open but also flat and closed 
FLRW solutions and one of the solutions provides the self-accelerating expansion of the late-time universe. 

Although the dRGT theory succeeds in providing the self-accelerating solution for the universe, the number of 
propagating degrees of freedom in the theory is not correct for the cosmological solution. Basically in four-dimensional 
spacetime, there are 5 degrees of freedom for the massive gravity while there are only 2 propagating degrees of freedom 
if the cosmological Friedman-Lemaitre-Robertson-Walker (FLRW) solution is assumed Moreover, De Felice, 

Gumrukcuoglu and Mukohyama 0 found that in the FLRW limit of the anisotropic solutions of the dRGT massive 
gravity, there is a non-BD ghost instability. This suggests the motivation for modifications on the dRGT massive 
gravity in order to cover the viable cosmological solutions. 

One of possibilities to obtain a viable cosmological model of massive gravity is that, one may generalize the dRGT 
theory by breaking the isotropy and study the anisotropic solution |15l . [ifij. An alternative possibility can be obtained 
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by introducing new degrees of freedom along with its coupling to the massive graviton. For example, one can promote 
the graviton mass to be a function of an extra scalar field, the so-called mass-varying massive gravity [T^ . One can 
also multiply a fiducial metric with the conformal factor depending on an extra scalar field |l 8 l |. The action of this 
extended dRGT massive gravity is invariant under quasidilaton global symmetry and then known as the quasidilaton 
dRGT massive gravity. Unfortunately, It was found that the self-accelerating solutions of the model are always 
plagued by ghost instability Thus, the quasidilaton massive gravity is further extended in order to avoid 

the ghost instability by introducing a co upling between the quasidilaton scalar field and Stiickelberg fields while the 
quasidilaton symmetry is still preserved |2ll - l23l | . It was also found that this extension can provide the correct number 
of degrees of freedom in the theory. It also was found that the effect of gravitational wave in quasidilaton massive 
gravity can be made larger than the Hubble parameter leading to an explanation of the suppression of the power 
spectrum in cosmic microwave background |24| . Recently, another possible way to obtain viable model of massive 
gravity is that the extension of dRGT massive gravity by introducing nontrivial coupling to the matter field [2^, [2^ . 

Another interesting extended massive gravity theory is Dirac-Born-Infeld (DBI) massive gravity [l^l . This extension 
is obtained by introducing DBI scalar field into the theory through the coupling term in such a way that the scalar 
possesses generalized Galileon shift symmetries . Like original dRGT massive gravity, it was shown that the model 
does not admit flat and closed FLRW solutions as well as the number of degrees of freedom is still not correct (^ . 
This may occur from using the Minkowski fiducial metric. Similarly to the original dRGT massive gravity, one may 
consider the generalized version of the fiducial metric in order to obtain the solutions which admit the flat and closed 
FLRW universe. Since the flat Minkowski fiducial metric in DBI massive gravity is motivated from the induced metric 
from brane world scenario, it is worthwhile to put a more general induced metric in consideration. In this work, we 
consider the extended dRGT massive gravity in which the fiducial metric is obtained from the induced metric of 
the brane moving in five-dimensional Schwarzschild-Anti-de-Sitter (Schwarzschild-AdS) spacetime [s^- It is found 
that the model not only admits the open but also flat and closed FLRW universe. The background solutions can 
be divided into two branches; a self-accelerating branch and a normal branch. For the self-accelerating branch, the 
graviton mass will play the role of the cosmological constant to drive the late-time expansion of the universe. For 
the normal branch, the general solutions are complicated and it is not easy to obtain the analytical solutions. To 
extract some information of this branch, we consider a specific class of the solutions and found that they can also 
provide the late-time expansion of the universe. The cosmological perturbations of the model are expanded around 
these background solutions. The results of cosmological perturbations in flat FLRW universe is that there are no 
ghost in the linear level while the number degree of freedom is still not correct in the self-accelerating branch. For 
the normal branch, there has the correct number of degrees of freedom but the ghost instability always presents in 
the theory, at least at small scale. This problem may be alleviated by introducing the coupling term similar to one in 
the generalized quasidilaton massive gravity. 

The paper is organized as follows: In Section El we present the setup of the model and find the background 
solutions. In Section uni the cosmological perturbations are analyzed in order to find whether the number of the 
degrees of freedom in the theory is correct and whether there is a ghost among them. Finally, we summarize the 
results in Section ITVl 


II. THE MODEL AND THE BACKGROUND EQUATIONS 


In this section, we will consider dRGT massive gravity including DBI Galileon term as follows 


S = SdRGT + •S'dbi, 

j\/f2 r 

=—^ J d,'^^V^[R[9]+‘^ml{C2[g,g] + a3£3[g,g]+a4C4[g,g])] 


where 5 ^ 1 / is the physical metric which defines the measurement on the spacetime and is the fiducial metric 
introduced to construct nontrivial nonlinear interaction terms in massive gravity. The first term in the action, the 
Einstein-Hilbert term corresponds to the Einstein’s general relativity. The second term is the interaction terms or 
mass terms, involving both the physical and the fiducial metric, characterized by 3 parameters; mg, aa, 04 , where mg 
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can be interpreted as a mass of graviton. The interaction terms are constructed, to avoid BD ghost, as follows Q, 


A[5,5] = i([/C]2-[/C2]), 

(2) 

>C3b,5] = ^([/Cp-3[/C][/C2]+2[/C3]), 

(3) 

^ ([/C]4 - 6[/C]2[/C2] + 3[/CT + 8[/C][/C3] - Q[1C^]) , 

( 4 ) 


where the square bracket denotes the trace operation with respect to the physical metric and the building-block tensor 
is defined as 


= ( 5 ) 

The square root denotes the tensor which upon being squared equals the Moreover, the hducial metric 

is an induced metric from the five-dimensional fiducial metric gAB, 

gp.u = dfj,X^d^X^gAB, ( 6 ) 

where all of X^ = the so-called Stiickelberg fields, transform as scalars introduced to restore 

general covariance of the theory. Lastly, the last and extension term is a leading order term of the action which is 
invariant under the Galileon shift transformation, also known as DBI action, where A is interpreted as a tension on 
the brane. 

For the physical metric, we will consider the Friedmann-Lemaitre-Robertson-Walker (FLRW) metric, 

+ a^{t)Q,ij{x^)dx'‘dx^. (7) 


Here, the latin indices run over all the spatial indices; 1, 2, and 3. The 3-dimensional tensor Oy is a 3-space metric 
defined by 




1 - 


( 8 ) 


Here, the curvature of the 3-space is defined by the value of k where the closed, flat, and open geometry correspond 
to K being positive, zero, and negative respectively. 

The dRGT massive gravity with flat Minkowski fiducial metric including DBI Galileon term have been investigated 
in It was shown that the model does not admit flat and closed FLRW solutions as well as the number of 

degrees of freedom is still not correct. These behaviors are similar to the dRGT massive gravity. In order to obtain the 
flat and closed FLRW solutions, one may generalize the fiducial metric from the Minkowski one. The generalizations 
in order to obtain the flat and closed FLRW solutions in dRGT massive gravity were investigated by using both FLRW 
and de Sitter fiducial metric [iMl- For the DBI massive gravity model, the fiducial metric can be interpreted as 
an induced metric in the brane world scenario. In order to generalize the fiducial metric to obtain the closed and 
flat solutions for the DBI massive gravity, we consider the form of the fiducial metric, gAB, which corresponds to a 
domain wall moving in five-dimensional Schwarzshild-AdS spacetime [s^, 


ds^ = -f{X^)dT^ + + j^^idX^f, (9) 

where 

= + ( 10 ) 

and K, g, as well as I are parameters of the model. The parameter k characterizes the three-dimensional surface 
corresponding to sphere for k = -|-1 for plane for k = 0 and hyperboloid for n = —1. The parameter I characterizes 
the curvature of the bulk spacetime and the parameter g characterizes the radius of a black hole in five-dimensional 
spacetime. Note that all parameters as well as X^ are dimensionless. In order to compare this fiducial metric to the 
physical one, it is convenient to rescale X^ to be (fiX where A is a parameter which has mass dimension. Then, the 
fiducial metric can be rewritten as 




( 11 ) 








4 


where 




fiX^) 

A 2 


K /i 

“ 3?^ 


£ 

P ■ 


( 12 ) 


Note that are Stiickelberg fields and (j) is an extra physical field. It is important to note that the additional scalar 
degree of freedom (j) is distinguish from the BD scalar since the BD scalar is always eliminated by construction of the 
mass terms. With general form of the fiducial metric, it was also shown that the BD ghost is not presented by using 
the Hamiltonian formulation even though the fiducial metric involves derivatives of the scalar field [s^. For the 
Lagrangian formulation, by usin g th e fiducial metric to be in the flat FLRW form, the scalar mode corresponding to 
the BD scalar is also eliminated [i3|. In this work, the elimination of a scalar mode corresponding to the BD scalar 
will be discussed in section UlI Cl 

In the brane point of view, this metric plays the role of the induced metric on the brane. By using this ansatz 
and unitary gauge, (p^ = t and = x®, the dRGT action in ([1]) can be expressed as 


SdRGT = iMh I d^xVUa^N ^+ml (f(X) - G(X)^)^ , 

«5'dbi = —SMpjTOgOfA J d^xV^a^NX^^, 

where D denotes the determinant of Q.ij{x^) and 

F{X) = {l-X){2-X) + ^{l- Xf (4 - X) + ^ (1 - X)" , 
G{X) = (1 - X) + a3 (1 - xf + ^ (1 - X)^ 

and we have defined 


(13a) 

(13b) 


(14a) 

(14b) 




a\ = 




(15) 


In order to obtain the equations of motion, one can vary the action in Eqs. (I13al) and (I13bl) with 
a. By varying the action with respect to N, we obtain 


respect to N and 


3-^Pi — Pg = —^MpiTjfF. 


(16) 


From this equation, one can see that the effective energy density, pg, contributed from graviton mass is dependent 
on time through the function F and it vanishes when rrf = 0 or X = 1. Varying the action with respect to a, the 
equation of motion can be written as 




N 


— Pg = —ZMpiwf 


F-^il-r] 


(17) 


where 


n = JL = JL- 

~ aN' ^~NX' 


(18) 


By considering the definition of pg in Ea. (ll 6 p and the effective pressure contributed from gravitational mass, Pg in 
Eg. (1171) . the solutions for accelerated expansion of the universe, Pg = —Pg, can be obtained by ^'(1 — r) = 0. We will 
see below that F' = 0 corresponds to the solutions in self-accelerating branch and r = 1 corresponds to a particular 
class of normal branch. 

For the equation of motion corresponding to variation of the action with respect to the Struckelberg fields, we will 
use the method investigated in (l^ by expanding the action up to the linear order in Sp^ = and 6cj) = tt® without 
variation of the physical metric. As a result, the first order perturbation of the action can be written as 


SS = ZMpirnf 


d'^xy/^a^N 




nF' {H - XH^)n° , 


(19) 
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where 



(j)n 


{G + aAX^) (6/3 + + / ( 2 # - 

27 ;? 


+ / - NnX 1 + 




A/ 



= 


jL 

(pn 


Therefore, the equations of motion can be written as 


J = 0, 

F' {H - XHP = 0. 


( 20 ) 

( 21 ) 


( 22 ) 


Now we have four equations of motion and we have three variables; a, N and </>. However, four of them are not 
mutually independent due to the Bianchi identity corresponding to the conservation equation. This constraint can be 
written as 


.5S ^df5S\ 

\Jn) 



_d n5S_ 

dt \n Sn 


= 0 . 


(23) 


From Eq. (1^^ . one can classify the solutions into two branches. The first branch corresponds to F' = 0 and the 
second branch corresponds to {H — XHp = 0. For the first branch, the equation can be expressed as 


F' = (-3 + 2X) - a3(l - X)(3 - X) - 0 : 4(1 - Xf = 0. 

Solving this equation, one obtains 

. , „ , 1 + 203 + 04 ± \/l + 03 + o| — 04 

X(j) = X±a, where X± = -- - -, 

0:3 + CH4 

which correspond to self-accelerating solutions since they provide 

Pg = —Pg = MpiA± 

= ((1 + a 3)(2 + 03 + 20 ^ - 304 ) ± 2(1 + 03 + o^ - 04 ) 3 /^) 

This branch of solutions are of the same expression with those from the dRGT massive gravity model with FLRW 
fiducial metric [l^- Note that the DBI scalar field has no contribution in this self-accelerating branch which should 
not be so surprising since such contribution to the gravity sector must be introduced via the minimal coupling between 
the scalar and the massive graviton like in the mass-varying massive gravity 0 or the quasi-dilaton massive gravity 
[ 1 ^ . Note that the solutions we have so far are valid for any kind of geometry; k can be set to any value while the 
setup studied in Ref. [ 2 ^ and in the pure dRGT theory [^, [l^ allows only the open slicing of the FLRW geometry. 

Considering the self-accelerating branch, one can substitute the self-accelerating condition, F' = 0, into Eq. ((22l) 
and then redefine variable as p = tp{(j)). As a result, Eq. (l22ll becomes 


(24) 

(25) 

(26) 


+ + + (27) 

The solution for this linear differential equation is 

0“=V- = £-(^)'/’, (28) 

where (pQ is an integration constant. From this solution, we can find the relation of the lapse function in terms of the 
DBI scalar field by using Eq. (HU) as follows 


iV2 


3A2 

{\^pK± - ?,kXI) 




(29) 
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Moreover, we can investigate the possible interval of the DBI scalar by considering (j)^ > 0. Its interval is determined 
by A, p, I, K and ^o- In case of the flat universe, it can be written in the simple form as 


0 < (j)< 


¥7^ 



(30) 


where b = From Eq. (1^51) . one can see that the hxed points of the system are / = 0 and / = The 

scalar field will evolve to the points = 0 or cjp = b for the fixed point / = 0. For the fixed point / = the 

/ _xl/2 

scalar field will evolve to the points (1 ± + . For the self-accelerating branch, it seems like 

that the dynamics of the universe is not controlled by the property of the fiducial metric since the Hubble parameter 
H is always constant and does not depend on function f. In other words, no matter how the scalar field tp evolves, 
the Hubble parameter H is always constant. 

For the second branch, namely normal branch, the equations of motion are very complicated. However, we can 
restrict our attention to a particular class of the solutions which simplihes the calculation while still provide the 
significant result of the system. The characteristic equation of this branch is H — XH^ = 0, so that we obtain 
X = H/H^. From Eq. (1^^ . we choose the solution such that G + = 0 to simplify the calculation as well as 

to capture the significant dynamics of the universe. This equation also satisfies the equation obtained by varying the 
action with respect to n, dS/Sn = 0, and can be expressed in terms of X as 

(a 4 — Scta) X^ — 3 (as -I- a 4 ) X^ -|- 3 (1 -I- 2 a 3 -I- 0 : 4 ) X — (3 -I- 3 a 3 -I- a 4 ) = 0. (31) 

Solving this equation, one obtains the relation between cp and a as 


X(j) = csa, (32) 

where C 3 is a constant depending on a 3 ,a 4 and a a. Since G + cxaX^ = 0 is the third order equation, C 3 can take 
three values to satisfy this equation. There is a particular value of aA such that 04 = 3q:a which provides only two 
solutions of X. Substituting this relation into the equation H — XH^ = 0, one obtains 


n 


N 


= C3- 


(33) 


This solution corresponds to r = 1. Note that all solutions where X is constant in normal branch provide the 
constraint of r = 1. As we have mentioned before, from Eqs. dm) and one can see that the constraint of r = 1 
provides the relation Pg = —pg = constant, corresponding to an accelerated expansion of the universe. This solution 
provides the same behavior like the solutions in self-accelerating branch. However, the solutions are different since the 
constant C 3 is not generally equal to X±. Moreover, the dynamics of the scalar held (p in both branches are different. 
In order to hnd the dynamics of the scalar held in normal branch for this solution, one can substitute relations in Eq. 
(|32j) and Eq. (|33j) into Eq. CH). As a result, the equation for the scalar held can be written as 


where 


Cr^cP^P 
1 -I- CnP"^ ’ 




Pg 

SMpicl 


(34) 

(35) 


Note that we considered this solution in the hat geometry where k = 0 for simplicity. This equation can be viewed as 
autonomous equation which has a hxed point at / = 0 corresponding to the point = b. To obtain the full behavior 
of the scalar held, one can solve Eq. (IMl) and the solution can be expressed as 


\/l + CrJ) tanh ^ 



\/l— Cnb tanh ^ 



25y^ Cn 


t + G, 


(36) 


where C is an integration constant determining the initial value of the held when initial time is taken. From this 
expression, we found that ^ b as t ^ 00 . This is what we expect from the analysis in the autonomous system. 
It is important to note that our analysis for normal branch is only a particular class of the solutions in the branch. 
Other solutions which may provide the result such that X = X[t) is much more complicated than this solution and 
thus is not easy to study analytically. They may give some interesting results and we leave this issue for further work. 
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III. PERTURBATIONS 


To investigate the stability of the model, we expand the action perturbatively up to the quadratic order. Here, the 
physical metric around which the action is expanded quadratically is expressed as follow. 


T ^9^^v^ 


(37) 


where the 


(0) 

9W 


is a flat FLRW metric and the perturbations are decomposed as 


5gw = 

Sgoi = N a{Bf + diB), 




+ d,Ej) + 25,,^ + 




(38a) 

(38b) 

(38c) 


where $, B, T, E are scalar parts, Bj, Ef are transverse vector parts, and hj^ is a transverse-traceless tensor part. 
We choose the unitary gauge in our analysis which corresponds to setting = 0. Thus the other perturbation comes 
from the DBI scalar field, S(j)- For simplicity, we choose to work in flat FLRW universe. In the following we will 
investigate the stabilities of the perturbations in tensor, vector and scalar modes separately. 


A. Tensor modes 


For the tensor modes, we expand action in Eq. © up to second order of . After that we keep only the second 
order and then transform the perturbation variables to those in the Fourier space. As a result, the second order of 
the action for the tensor modes in Fourier space can be written as 


where 


S{ 2 )^Mk I Skdta^N 



Mqw — 


A = X 





+014 (1 



(2-A)r^ 


(39) 


(40a) 


(40b) 


The kinetic term of the tensor mode always has the correct sign. Thus there is no ghost in this mode. The absent of 
tachyonic instability is required by the condition > 0. For the self-accelerating branch, this condition 

depends only on 0:3 and a 4 . We use this condition together with pg > 0 and X± > 0 in order to find compatible 
regions in the ( 03 , 014 ) space. The allowed regions for both solutions are shown in Fig. [T] The left panel corresponds 
to X^ solution and the right panel corresponds to A_ solution. The shaded region with horizontal-dashed-blue line 
correspond to the condition > 0- The shaded region with vertical-dashed-black line corresponds to the condition 

X± > 0. The grey region corresponds to the condition pg > 0. For the condition > 0, it is found that the 

region will exist if r > 1 for X+ solution and r < 1 for A_ solution. For r = 1, it corresponds to = 0 which is 

reduced to massless gravity theory. Since the solutions in this branch are n ot g overned by the DBI-scalar field, the 
region are the same with dRGT massive gravity with FLRW fiducial metric [1^ . 

For the normal branch, we choose the particular class of the solutions which discuss in the previous section. These 
solutions depend on three parameters, 013 , 014 and a a. For simplicity, we also restrict our attention only in the case 
ttA = 014 / 3 . This restriction provides two solutions by solving Eq. (I^Tl) . As a result, the allowed region which satisfies 
the condition Mq-^ > 0, pg > 0 and A_ > 0 is shown in Eig. Note that there is no allowed region to satisfy the 
condition pg > 0 for the solution A+. Surprisingly, the graviton mass for these solutions coincides with the effective 
mass investigated by using de Sitter and FLRW fiducial metric in dRGT massive gravity 0 - This effective mass 
is constrained by Higuchi bound as Mq^ > 2H^. We use this condition together with pg > 0 and A_ > 0 to find 
the allowed region. The result is shown in Fig I2bl At this point, one can see that we can find the viable model 
by specifying the parameters as and a 4 in both branches. However, by investigating the perturbation in vector and 
scalar mode, it is found that the theory is suffered from ghost instability. We explore this behavior in the next two 
subsections. 







as 


as 


(a) X+ solution 


(b) X- solution 


FIG. 1: Allowed regions in space for the self-accelerating branch plotted by using three conditions: 

Mqw > 0, Pg > 0 and X± > 0. The left panel corresponds to X+ solution and the right panel corresponds to A_ 
solution. The shaded region with horizontal-dashed-blue line correspond to the condition Mq^ > 0. The shaded 
region with vertical-dashed-black line corresponds to the condition X± > 0. The grey region corresponds to the 

condition pg > 0. 



as 

(a) X- solution with condition > 0 


as 

Jlyf2 

(b) X- solution with condition > 2 


FIG. 2: Allowed region in ( 0 : 3 , 04 ) space for normal branch plotted by using three conditions, Mq^/ > 0, pg > 0 and 
X- > 0 for the left panel and Mq^/H^ > 2, pg > 0 and X- > 0 for the right panel. The shaded region with 
horizontal-dashed-blue line correspond to the condition Mq^ > 0 in the left panel and Mq^/H'^ > 2 for the right 
panel. The shaded region with vertical-dashed-black line corresponds to the condition X± > 0. The grey region 

corresponds to the condition pg > 0. 
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B. Vector modes 


For vector modes, we follow the calculation step as performed in tensor modes. As a result, the second order action 
in Fonrier space can be written as 


5'(2) = J dPkdta^N ■ 


\El?- 


8V2 


2Na 


-bTeT - - 


2H 

w 


\E 


■T\2 


2a2 


3H‘ 


bT?- — 


+ -M^ 

Ml, 2 ' 


GW 


eT?- 


Pg 

Ki 




( 1 - 


(41) 


From this action, one can see that Bf is non-dynamical. Thus we can algebraically solve it and then substitute its 
solution back to the action. As a result, the action becomes 


5(2) = ^ I dPkdta^N 


k’^nri^a^XF' 

(A:2(i + r) - 2mla^XF') 


\ET\^-YMgw 



(42) 


It is found that, for self-accelerating branch; F' = 0, there is no propagating d.o.f in vector modes. For the normal 
branch, it is ghost free if rn^F' < 0. In the case of m^F' > 0, the propagating vector seems to suffer from ghost at 
small scale or large k when 


fc2 > 


2a^XF' 

(1 -hr) 


■m" 


m 


9' 


(43) 


The sound speed of propagation can be written as 


Mgwi^ + r) ^ A{l + r) 
2mlXF' 2XF' 


(44) 


Thus, the condition to avoid gradient instability is 


A 

J' 


< 0 . 


(45) 


It can be summarized that the condition for avoiding the instabilities is that rUgA > 0, < 0 or corresponding to 

rn^gF' < 0. 


C. Scalar modes 


For scalar modes, the second order perturbations in Fourier space for Einstein-Hilbert can be written as 


S^eI = ^ f d^kdta^Nr^ 


F^ 




y 


18a2 

-h ik'^H^B^ - 9 i 72$2 


2H 


6^2 

yr 

-h3i?2^ 

2k^-^E 

3a2 


2k^BE 

SaN 


■ /12i/4> Ak'^B 
4' 


V N 


4-2 


4> 


v ® 

Uk^HB 

V a 


2k^ ^ 2H 

—2“ + 3 I 


4- 


N 

4fc2 


aN 


3iJ2 


18iJ^ 


2k^E 

3a2 


(46a) 
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For dRGT and DBI parts, the action can be expressed as 


c(2) C 

^dRGT^^ 


)--WTT7rj 




a?NnX 


+ 3 


Pn 


Mhml 


+ 2A\'^‘^ + k^ 


H - — 

n^F' 

MpiirP N‘^X{l + r) 


y4\ 2 

3 V Mpim^ ' 2 


Pn 




Pa _$2 


Mhml 




Nnp 

+ 18 ( F" - 1 $'1' 


iVA/(l + r) 


(46b) 


From the action in (|46al) and (I46b|) . B and 4) are non-dynamical. We can use their equations of motion to eliminate 
them from the action. These constraints comes from the energy and scalar part of the momentum conservation. Now 
we have three variables 5(t),E and 4>. One of them is a scalar degree of freedom for the massive graviton, the others 
are BD scalar and DBI scalar. However, we found that the kinetic term of 'k will vanish by imposing the background 
equation (fTO)) and dig. Therefore, the scalar 'k can be interpreted as the BD scalar since it must be eliminated by 
construction of the theory. We can eliminate it from the action by following the previous procedures. Note that the 
full lagrangian before integrating out 'k is quite lengthy and then we left in the Appendix Ea. dAll) . Generally, we 
have the scalar mode action of the form 

S = j d^kdt {Kijx^x^ + Mijx^x^ + PiJx'x^) , (47) 

where xi = ^4' s-nd X 2 = E- The Kij,Mij and Pjj each are elements of corresponding 2x2 matrices associated to 
each combinations of the remaining propagating d.o.f. 

For the self-accelerating branch, we have F' = 0 = H. It turns out that according to this branch, there is no 
propagating scalar d.o.f in gravity sector; only propagating scalar is obviously the DBI scalar. The scalar mode action 
for the self-accelerating branch is specified with the matrix Kjj as follows. 


_ iMpiiTiga^ [G + a\X^) 
" 2n3 


det(Ff|S)) = 0. 


which the only nonvanishing matrix element is Kn corresponding to the DBI scalar Sif) as claimed. 
For the normal branch, we consider the case which satisfies the condition we considered previously. 


G + aAX^ = Q, r=l. 


(48) 


(49) 


Even though this equation is not independent to the other equations of motion, we still use this equation in order to 
simplify the second order action. By imposing this background equation of motion, it seems like that the DBI-scalar 
d.o.f does not propagate. However, substituting back the solution of 41, one obtains 


= - 


QMhmlaX^V'^F''^ 

mpEw ' 


-^12 — 


k: 


(n) 


3Ml.ik^XVF'^ 
NH^UW ’ 
k^a^F' 


— ^pi 


k^F'^ 


mu miaNH^U^W) 


detiK\y) = - 


QMl^imlk^a^X^V^F'^ 

mrnnJw ^ 


(50) 


where 


U 

w 


3F'- 


4fc2 


m^a^X 


V= 1 


12A-k2 


2fci 






, , 3mlxm 
+ 9X]F' + - ^ - 


iF2 



(51) 
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It is found that, for the particular class of the solutions in normal branch, the d.o.f of the theory is still correct. At 
small scale, U oc and then W oc —F' where the proportional constant is positive. Substituting these results 

into it is found that oc m^F' where the proportional constant is positive. This leads to the condition 

to avoid ghost such that rri^F' > 0. This condition is in conflict with the condition from the vector mode in which 
rn^gF' <0. As a result, the ghost degree of freedom always exist in this branch at small scale (either in vector modes 
or in one of scalar modes). For the other mode of scalar perturbations, one can find the result by substituting U and 
W at small scale into det . As a result, we found that det is always negative implying that, at least, one 

of scalar mode is always ghost. It is important to note that the leading term of is proportional to k~'^ leading 

(n) 

to —>■ 0 when we take k oo. This behavior occurs only in the particular solutions we choose to consider 

since the leading term, corresponding to is proportional to G + ■ This suggests that the contributions from 

linear perturbations are very small. Therefore non-linear perturbations may have nontrivial contributions and the 
conclusion of our results may change when non-linear perturbations are taken into account. 


IV. CONCLUSION 

In this work, we have considered the extended DBI massive gravity by generalizing the fiducial metric. In the 
original version of DBI massive gravity, the fiducial metric is in Minkowski form. It was found that the flat and closed 
solutions for FLRW metric are not admitted. Moreover, the number of degrees of freedom of the massive graviton is 
not correct. The flat Minkowski fiducial metric in DBI massive gravity plays the role of the induced metric on the 
brane. In the present work, we generalize the fiducial metric by considering the induced metric corresponding to a 
domain wall moving in five-dimensional Schwarzschild-AdS spacetime. 

From background equations, we found that the solutions not only admit open geometry of FLRW metric but also 
flat and closed geometry. Moreover, the solutions can be divided into two branches; self-accelerating branch and 
normal branch. For self-accelerating branch, the graviton mass plays the role of cosmological constant to drive the 
late-time acceleration of the universe. We found that the Hubble parameter H is always constant and does not 
depend on function / which characterizes the form of the fiducial metric. Therefore, the dynamics of the universe 
are not controlled by the property of the fiducial metric. Furthermore, the DBI scalar field has no contribution 
to the solutions in this branch. This is not so surprising since such contribution to the gravity sector must be 
introduced by the minimal coupling between the scalar and the massive graviton like in the mass-varying massive 
gravity and quasi-dilaton massive gravity. In order to investigate the number of degrees of freedom of the theory, we 
performed cosmological perturbations around this background solutions. From the tensor modes, we have explored 
the conditions to avoid the tachyonic instability, Mqt^ > 0, while the ghost and gradient instabilities are absent. By 
using this condition together with consistent conditions, Pg > 0 and X± > 0, we found that there are allowed regions 
in (as, a 4 ) space in both A+ and A_ solutions as shown in Fig. [T] By including the investigation of vector and scalar 
modes, we found that the number of degrees of of freedom is not correct. There are only two degrees of freedom 
for gravity sector contributed from tensor modes while the contribution from vector and scalar modes vanish. This 
inconsistency is similar to the dRGT massive gravity and DBI massive gravity. 

For the solutions in normal branch, the equations of motion are complicated and it is not easy to obtain the 
analytical solutions. Therefore, the solutions in this branch are quite lack of investigation in literature. This is not 
only due to complicated equations but also the solutions generally do not provide an accelerated expansion of the 
universe. In the present work, we restrict our attention to a particular class of the solution in order to simplify the 
calculation and also to capture some significant behaviors of the solutions. A characteristic constraint of this class of 
the solutions is r = n/{NX) = 1. We found that this class of the solutions also provides the accelerating expansion of 
the universe like in self-accelerating branch inferred from the definition of effective energy density pg and the effective 
pressure Pg in Eqs. m and (fT7l) respectively. Even though the solutions in this class provide the same accelerated 
expansion of the universe like the solutions in self-accelerating branch, these solutions are different since the solutions 
in this branch depend on the effect of DBI scalar field via the coupling cta. By setting the coupling = a^/i to 
simplify the solutions, the allowed regions plotted by using conditions pg > 0, A_ > 0 and Mqt^ > 0 are explored. 
Moreover, we also found that the graviton mass in this branch coincides with the effective mass investigated by using 
de Sitter and ELRW fiducial metric in dRGT massive gravity 0. This effective mass is constrained by Higuchi 
bound as Mqt^ > 2H^. We use this condition together with Pg > 0 and X- > 0 to find the allowed region in (as, a 4 ) 
space shown in Eig [21 By using the investigation of the perturbation in vector and scalar modes, we found that 
the number of degrees of freedom is correct. However, at small scale with > rri^, the conditions to avoid ghost 
instability from vector modes and scalar modes are in conflict implying that there exists at least one ghost degree of 
freedom. This investigation suggests that the extended model of DBI massive gravity in this way dose not provide a 
viable cosmological model to explain the late-time acceleration of the universe. In order to obtain a viable extended 
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model of DBI massive gravity, one may introduce the coupling term similar to one in the generalized quasidilaton 
massive gravity or introduce the nontrivial coupling matter helds like in original dRGT massive gravity. 
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Appendix A: Pull quadratic scalar mode action 

The full quadratic scalar mode action, after imposing the background equations of motion (HSlini), is 


S'(2) = Mh I d^kdta^Ni 


3to2 f G + uaX^ 
Nn^ 


Scj)^ 


k^F' 




k^rF' 

a?NHU 






aNH 




18 \ 2a2 4 ' 2a^mU 

^ /c2(G' + 3aAA2) 3J, 


a^NXfXHU 
k^r'^p’ k^ F r)(j)U F + r)(j) - (j)r^U^ 


mla'^Ncj^XHU'^ 


E^ 


m 


F' (2k'^{f F N'^rX'^) 9mlN^X^F''^ 


a^N^rX^ 


a^N 20 ?^ \ /r(l -f r)A3 


mu 


+ 


6fc2iV2<^riJ (2A2/(1 -h r) -f a^rXHH^) F '' 
X^pH{lFr)U ^ 

k^{2PrHH^ F imlX^fF')F' 


„ mj,W „ 
5p -f ^—4-^ 


Qa^XfmU 


kP(2P F mlPXU)F' 

6Pmu 


'i/E 




GW 


+ 


F' 


inPgXX 


1 - 


' 6m2A 2kUrH 3mlf'{p F N^pX^Hl)' 

a a^XfHU 2mprX j 

XrHA 3m2 (aV(^ + ^^(1 - r)U) - 2PXrHH^) F'‘- 


H 


F" - 


2aXfmU 




(5(^4' 


(Al) 


2amU 

By solving 'll and then substituting back into the action, the components of the kinetic terms matrix can be written 


All = - 


9MPm^aX^V^F'^ SAIpiiriga^ (G -I- 


?,Mlik‘^XrVF''^ 


Ko2 = XP 


xmw 

PPF' 


2n^ 


PPF'^ \ 


I 4AG rnlaNmmW j 


det A/j = — - 


= 

xmuw 

9A/|,;m2fc4a4y2y2^/3 

4mmuw ’ 


(A2) 

(A3) 


where 


U = 

W = 


3F' - 


2P (iFr) 


rn^a^ 


12A-f 2 


A 

2PP 


rrigpH'^U 


V = [IF 


6A 


aH(j)HA 

A/ J’ 


3 XPH 4 

H 


F' F 


3mJ,X^rF'^ 


m 


- 6AM 1 




(A4) 
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